Bubble-bubble interaction: A potential source of cavitation noise 
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The interaction between microbubbles through pressure pulses has been studied to show that it 
can be a source of cavitation noise. A recent report demonstrated that the acoustic noise generated 
by a shrimp originates from the collapse of a cavitation bubble produced when the shrimp closes 
its snapper claw. The recorded acoustic signal contains a broadband noise that consists of positive 
and negative pulses, but a theoretical model for single bubbles fails to reproduce the negative ones. 
Using a nonlinear multibubble model we have shown here that the negative pulses can be explained 
by considering the interaction of microbubbles formed after the cavitation bubble has collapsed 
and fragmented: Positive pulses produced at the collapse of the microbubbles hit and impulsively 
compress neighboring microbubbles to generate reflected pulses whose amplitudes are negative. 
Discussing the details of the noise generation process, we have found that no negative pulses are 
generated if the internal pressure of the reflecting bubble is very high when hit by a positive pulse. 

PACS numbers: 47.55.dp; 47.55.dd 
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I. INTRODUCTION 

In a recent paper, Versluis et al. reported that the snap- 
ping shrimp {Alpheus heterochaelis) living in the ocean 
can generate a cavitation bubble by rapidly closing its 
large snapper claw The rapid closure produces a 
negative pressure in seawater, by which cavitation nu- 
clei (e.g., air microbubbles) are explosively expanded to 
a radius of a few millimeters. The cavitation bubble 
then collapses violently and emits a loud acoustic noise. 
The experimentally recorded acoustic signal presented in 
Ref. 1| consists of a strong (positive) pressure pulse, 
clearly produced at the bubble collapse, and a subse- 
quent broadband noise A single-bubble theoretical 
model (the Keller-Miksis equation) succeeded to repro- 
duce the strong pressure pulse (and a weak precursor 
signal) but failed to describe the broadband noise. Ver- 
sluis et al. stated that the broadband noise is, partly, due 
to the reflection of the pressure pulse at nearby aquarium 
walls. However, the broadband noise begins earlier than 
the reflected wave reaches, immediately after the bubble 
collapse. 

The broadband noise appears to consist of both posi- 
tive and negative pressure pulses (or steep spikes) whose 
amplitudes are smaller than the first strong pulse. As 
demonstrated in Ref. [IJ and shown below, however, the 
single-bubble model cannot describe negative pulses. A 
key to resolving this inconsistency is given from an image 
recording of the cavitation bubble. In a series of high- 
speed images, it was found that at collapse the single 
cavitation bubble breaks apart through the surface insta- 



bility and then an opaque cloud of microbubbles appears 
The bubble cloud seems to grow and finally dissolve 
away. 

We hypothesize that the interaction between the mi- 
crobubbles through pressure pulses is a source of the 
negative pulses involved in the broadband noise. It is 
known that bubbles undergoing volume change interact 
with each other through the pressure waves that they 
emit. Bubble-bubble interaction of this type leads to a 
variety of phenomena, including attraction/repulsion be- 
tween bubbles [1, |3, d, 0, H, Is], filamentary structure 
formation [l(]||. change of eigenfrequencies [ll|, superres- 
onances (l2j|. emergence of transition frequencies [§, [l^ , 
avoided crossing of resonance frequencies sound lo- 
calization and suppression of cavitation inception 
(TgI . [13] ■ Because at generation the microbubbles are 
highly compressed, they must begin volume change (as 
observed) and emit pressure pulses at their collapse. Tak- 
ing the bubble-bubble interaction through the pressure 
pulses into consideration, in this paper we suggest a pos- 
sible origin of the negative pulses. Here we do not aim 
at providing a quantitative explanation because the ac- 
tual sizes and number density of the microbubbles are 
now unknown. We instead attempt to elucidate the ba- 
sic mechanism of negative pulse generation, reducing the 
problem to the interaction of only two bubbles. 

The rest of this paper is organized as follows: In Sec.lTTl 
the model equations and assumptions used in this study 
are introduced. Section IIIII presents numerical and the- 
oretical results and discussions on the pressure pulses 
emitted by interacting bubbles, and Sec. IIVI summarizes 
the obtained results. 
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II. MODEL EQUATIONS 

The theoretical model used in this study is the coupled 
Keller-Miksis equations [l3| , which describe the radial 
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motion of two coupled spherical bubbles in a compressible 
liquid: 
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where ri is the distance measured from the center of bub- 
ble i. This pressure equation can be given from the 
continuit y a nd Euler equations of fluid flow (see, e.g., 
Refs. and is in the following used to examine 

the acoustic signal from the bubbles. Time-delay effects 
[23I I due to the finite sound speed of water are neglected, 
but a remark will be made on a consequence of it. As a 
first approximation we neglect the translational motion 
of bubbles due to the secondary Bjerknes force (an inter- 
action force proportional to that is, assuming 
dDij/dt = 0. Hence, we set D^- to be much larger than 
RiQ {Dij « 10i?io) and consider only the first few periods 
of bubble oscillation, in which the translational velocity 
should be very small. As demonstrated below, the first 
one or two periods are sufficient for our discussion. 



where Ri — Ri{t) is the instantaneous radius of bub- 
ble i, Dij is the center-to-center distance between bub- 
bles i and j, and the overdots denote the time derivative 
d/dt. The surrounding liquid is assumed to be water 
with density p = 1000 kg/m^, viscosity ^ = 1.002 x 10^^ 
kg/ (m s), sound speed c = 1500 m/s, and surface tension 
a ~ 0.0728 N/m. The far-field pressure Pq is assumed to 
be constant in time and equal to the atmospheric pres- 
sure, 0.1013 MPa. The gas in the bubbles is assumed 
to be a van der Waals gas (air) and the pressure inside 
bubble i {pb,i) is determined by 
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where Rio is the ambient radius and hi is the hard-core 
radius (i?io/8.54 for air [11]). The polytropic exponent 
of the g assumed to be equal to its specific heat 

ratio 7 (1.4 for air) because our interest is in the pres- 
sure pulses emitted at bubble collapse, where the bub- 
ble behavior is nearly adiabatic. We confirmed numeri- 
cally that even when the heat exchange between the bub- 
bles and water is taken into account by, e.g., a switch- 
ing method for k [l^, results are essentially the same 
as those shown below. The vapor pressure is neglected 
for simplicity. We do not consider mass exchange (i.e., 
evaporation/condensation, mass diffusion) and chemical 
reactions, which may occur inside the bubbles [20| . since 
they are not essential for sound emission: As demon- 
strated in earlier work on single bubbles [2l|, , models 
that do not take mass exchange and chemical reactions 
into account can describe bubble emitted pressure waves 
with sufficient accuracy. 

The last term of Eq. ((!]) describes the bubble-bubble 
interaction through the bubble emitted pressure waves 
and acts as a driving force on bubble i. This term was 
derived from the following simple formula, which corre- 
sponds to an equation for the pressure wave emitted by 
a pulsating sphere. 
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III. RESULTS AND DISCUSSIONS 

First we consider a single-bubble case to confirm that 
single bubbles can emit positive pulses only. The radius- 
time curves for two different bubbles of i?io = 10 pm, 
R20 — 15 pm with D12 = 00 are shown in Fig.[TJa). Here, 
we set the initial condition as Ri{t = 0) = 0.13i?io and 
Ri{t = 0) = in order to simulate the highly compressed 
state of the microbubbles at their generation. Due to the 
high internal pressure, the bubbles undergo rapid expan- 
sion and reach a maximum size of about 2.5i?io- Then 
they collapse and rebound many times. The oscillation 
amplitudes decrease monotonically in time due to the vis- 
cosity and compressibility of water. The pressure waves 
from the bubbles measured at ri ~ r2 = 30i?io are shown 
in Fig. [TJb). From this, one knows that the single bub- 
bles can emit positive pulses only. Between the positive 
pulses, one finds low-amplitude negative pressures, which 
resulted from a weak deceleration of the bubble surface 
acting when Ri{t) > RiQ. Their amplitudes (and also 
fundamental frequencies) are, however, obviously much 
lower than those of the positive pulses and hence they are 
incapable of explaining the negative pulses in the acous- 
tic signal. The positive pulses are produced at the bubble 
collapse where the bubble surface is strongly accelerated 
[i.e., Ri in Eq. ^ has a very large value]. However, from 
Eq. ^ one finds that to produce a negative pulse the 
bubble surface needs to be strongly decelerated: this may 
be impossible when only single bubbles are considered. 

This difficulty is resolved by considering bubble-bubble 
interaction. As well known, when a pressure wave prop- 
agating in water hits the water-air interface, most of 
its energy is reflected and a reflected wave is produced 
whose phase is opposite to the incident wave. This is 
because air's acoustic impedance is much smaller than 
water's. It is also known that negative pressure waves 
generated when strong positive pulses from collapsing 
bubbles hit a water-air interface can be strong enough 
to cause secondary cavitation [2^. Since the surface 
of gas bubbles considered here is also a free surface, it 
would be able to produce negative pulses. In Fig. ^ we 
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FIG. 1: (Color online) Bubble radii -Ri,2 (a) and emitted pres- 
sures pi, 2 measured at ri,2 ~ 30-Rio (b) in a single-bubble case 
{Di2 — oo) as functions of time. The ambient radii of the 
bubbles are Rio — 10 /xm and R20 ~ 15 /xm. The solid and 
dashed lines are for bubbles 1 and 2, respectively. The pres- 
sures are normalized by the atmospheric pressure Pq. Only 
positive pulses are found in pi,2. 



show a result for a double-bubble case. Here, we used 
the same parameters as in the above example except 
for D12 — 10i?io. In this case, the change of oscilla- 
tion amplitude is non-monotonic because of the modula- 
tion effect due to bubble-bubble interaction. In the bub- 
ble emitted pressures presented in Fig. OJb), in contrast 
to the single-bubble case, not only positive pulses but 
also negative pulses appear. In order to confirm the ro- 
bustness of this observation, we performed a parametric 
study. The result is presented in Fig. [31 where the max- 
imum amplitudes of the negative pulses for Riq = 10 
/ini and D12 — 10i?io are shown as functions of i?20- 
This proves that strong negative pulses are emitted in 
most cases. From this figure it is also found that the 
maximum amplitudes have a complicated dependence on 
the ambient radius. The radius-time and pressure-time 
curves for i?2o — 6 /im and R20 = 10 /xm are shown in 
Figs, m and respectively. For i?20 = 6 /xm, negative 
pulses are found but their amplitudes are very low. For 
i?20 = 10 ^m (i.e., i?20 = -Rio), no negative pulses are 
found, implying that systems of identical bubbles do not 
emit negative pulses. 

The frequency spectra for _Rio = 10 /im and _R20 — 15 
/xm taken from pi/Po and P2/P0 in Figs. [TJb) and[2l^b) 
are shown in Fig. [HI All of the presented spectra show 
a broad distribution. The peak frequencies in the single- 
bubble case are 0.287 MHz for pi and 0.185 MHz for p2, 
which are roughly the same as the repetition rates of the 
positive pulses deduced from Fig. [Hb) (about 0.27 MHz 



FIG. 2: (Color online) Same as Fig.[T] but in a double-bubble 
case (_Di2 = 10-Rio). Both positive and negative pulses are 
found in pi.2, as in the broadband noise reported in Ref. [l|]. 
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FIG. 3: (Color online) Maximum amplitudes of the negative 
pulses normalized by Po [min(pi,2)/fo] as functions of R2o- 
Rio and D12 are fixed to 10 fim and lOiiio, respectively. The 
circles and squares are for bubbles 1 and 2, respectively. 



and 0.17 MHz, respectively). We found several quali- 
tative differences between the spectra in the single- and 
double-bubble cases. One of the differences is clearly seen 
in Fig. [7l which shows the spectra in a frequency range 
around the peak frequencies. The spectrum of pi in the 
double-bubble case has a more complex structure than 
that in the single-bubble case, particularly in a frequency 
range between 0.1 MHz and 1 MHz: A number of char- 
acteristic peaks are added and the spectrum structure 
becomes much denser through the bubble-bubble inter- 
action. Compared to pi , the change in the spectrum of P2 
is small: No qualitative differences are found between p2 
in the single- and double-bubble cases. This suggests that 
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FIG. 4: (Color online) Same as Fig. (2] but for _Rio = 10 fim 



o 




Frequency [Hz] 




FIG. 5; (Color online) Same as Fig.El but for _Rio = R20 = 10 
fj.m. 



FIG. 6: (Color online) Frequency spectra of pi/Po (a) and 
P2/P0 (b) for 7?io = 10 ^m and R20 = 15 /im. The dashed and 
solid curves are for the single- and double-bubble cases, re- 
spectively. The spectra were generated using the data shown 
in Figs, mb) and[2l[b) in a time period from t = 2.5 \vs, to 100 
/is. The sampling frequency assumed is 500 MHz. 
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the larger bubble (bubble 2) has a more significant influ- 
ence on the neighboring bubble and the dynamics of the 
neighboring smaller bubble (bubble 1) is thus changed 
more drastically (the same tendency is found also in sys- 
tems of linearly oscillating bubbles in a sinusoidal sound 
field 0). When i?io = i?20 = 10 /im, the difference be- 
tween the spectra in the single- and double-bubble cases 
is very small, as can be expected from the above observa- 
tion that only positive pulses are emitted for i?io = i?2o; 
see Fig. |S1 The peak frequency is decreased by bubble- 
bubble interaction, which appears to be consistent with 



FIG. 7: (Color online) Frequency spectra of p\/Po (a) and 
P2/-P0 (b) for Tiio = 10 /xm and J?20 = 15 /im in a frequency 
range around the peak frequencies. 



the fact that the natural frequency of identical bubbles 
oscillating in phase each other is lower than that of the 
individual bubbles [HI, [S Ull ■ 

Let us consider how the negative pulses were produced. 
As can be seen in Fig. Hfb), the negative pulses from a 
bubble coincide with the positive pulses emitted by the 
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FIG. 8: (Color orrline) Same as Fig.El but for _Rio = -R20 = 10 
fim. 



other bubble at its collapse. This observation tells us that 
the negative pulses are produced when the positive pulses 
hit the surface of the neighboring bubble. Shown in Fig.lH] 
is a close-up view of pi and Ri in Fig. [2] around the time 
(t « 7 fis) of the first negative pulse from bubble 1. The 
figure also shows the inertia and acceleration portions of 
Pi, that is, 



and 



ai 
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One can see that the surface of bubble 1 is strongly de- 
celerated in the period when it emits the negative pulse, 
though no noticeable disturbance is seen in Ri because 
the duration of the deceleration is very short. This strong 
deceleration is clearly caused by the strong positive pulse 
from bubble 2, which impulsively compresses bubble 1. 
Deceleration is also observed in a following period where 
the bubble is shrinking, but no negative pressure is found 
because a large inertia (ai) cancels the effect from the 
deceleration 

In order to more deeply understand the mechanism of 
negative pulse generation, we perform theoretical studies 
based on the coupled Keller-Miksis equations. First, we 
consider the relation between the amplitudes of the inci- 
dent and reflected pulses. If Ri > Rio (i.e., pb^i < Pq) 
and i?i <C c at the time of collapse of bubble 2, where 
i?2 has a very large value, only the acceleration terms in 
Eq. ID) may be significant for bubble 1 and the remaining 
terms are thus negligible. Hence, the equation and also 
Eq. (lU are reduced to 
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FIG. 9: (Color online) (a) Close-up view of pi and Ri in Fig. [2] 
around the time of the first negative pulse from bubble 1, and 
(b) inertia {ot\) and acceleration (/3i) portions of pi in the 
same period, pi, qi, and /3i are normalized by Po- 



In Eq. ([7]), Ar is an acoustic Reynolds number defined 

here as 



A„ = 
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Since Ar is much larger than unity {Ar ~ 973 for Ri — 
2.6 /im, a typical minimum radius at collapse), Eq. ([7]) is 
further reduced to 
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If one measures pi and p2 at the same distance from 
the corresponding bubbles (i.e., ri = r2), Eq. pH)) is 
rewritten using Eq. ([8|) as 



Pi 
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Equation pO|) proves that Ri has the opposite sign from 
that of R2, that is, the surface of bubble 1 is strongly de- 
celerated at the moment when bubble 2 collapses, where 
the surface of bubble 2 is strongly accelerated due to its 
high internal pressure. From Eq. (jlip. one finds that 
the ratio between pi and p2 is simply determined by 
—R1/D12, and that the amplitude of the negative re- 
flected pulse is large if Ri is large when the pressure 
pulse from bubble 2 hits bubble 1. The numerical result 
shown in Fig.[2]^b) gives pi/p2 — -0.200 and -R1/D12 = 
—0.192 for the first negative pulse from bubble 1 (at 
t « 7 fis), and P2/P1 = -0.243 and -R2/D12 = -0.231 
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for the second negative pulse from bubble 2 {t !v 8 fis), 
both of which reveals a reasonable agreement. However, 
Eq. ([TT|) gives a less accurate result for the first negative 
pulse from bubble 2 (t w 5 ^s): P2/P1 — —0.507 but 
—R2/D12 — —0.369. This may be because pi and p2 are 
not large enough to fully satisfy the assumptions used in 
deriving Eq. PU)) . 

Next, we consider what occurs when i?io = i?20- As 
shown in Fig.[5l no negative pulses are emitted for i?io — 
R20 = 10 /im, although strong positive pulses are emitted 
which definitely hit the neighboring bubble. Here we 
attempt to explain this observation. For Rio = R20 and 
i?i = i?2, at the final stage of bubble collapse where 
Ri and pbA are very large and Ri ~ 0, Eq. ([T]) may be 
reduced to 



pRiRi 
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(12) 



where we neglected the surface-tension and viscous forces 
and Pq since their magnitudes should be much smaller 
than the internal pressure Since I/Ar (~ 1/977 for 
Ri — 2.61 fim at the first collapse) is much smaller than 
R1/D12 (~ 1/38), Eq. nil) is further reduced to 
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FIG. 10: (Color online) Noise generation process. Positive 
pressure pulses are emitted by collapsed bubbles (the right 
circle) and propagate through the surrounding liquid. The 
positive pulses generate negative reflected pulses when they 
hit and impulsively compress neighboring bubbles (the left 
circle). The positive and negative pulses will make up a cav- 
itation noise. Note that the pulse profiles shown represent 
those measured at a fixed location, as functions of time. The 
wavelengths of the pulses are much larger than the bubble 
radii. 



pRiRi — Pb,i — -^—R\Ri. 
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IV. SUMMARY 



From this we have 



Pb,i 



(14) 



This says that Ri is positive at the bubble collapse and 
hence the bubbles emit positive pulses. This conclusion 
is consistent with the above numerical finding and is not 
altered even if l/An is considered. 

This result comes from the fact that, in the considered 
case where the bubbles collapse at the same time, the 
sound pressure -jy—R\Ri from the neighboring bubble 
cannot exceed the bubbles' internal pressure Pb,i- The 
sound pressure generated by bubble 2, whose internal 
pressure is Pb,i^ must be smaller than when it is 
measured at the position of bubble 1, because the bub- 
bles are separated by a finite distance I?i2. Since the 
internal pressure of bubble 1 is also Pb,i, it is not ex- 
ceeded by the sound pressure, and the right-hand side of 
Eq. (|13p is thus positive: in other words, bubble 1 emits 
a strong positive pulse, resulting from the high internal 
pressure, whose absolute amplitude is greater than that 
of the negative reflected pulse. This finding suggests that 
the bubble surface, a free surface, docs not always pro- 
duce negative reflected pulses and the characteristics of 
the reflected waves depend on the state of the gas in the 
bubbles. The negative pulses found in Figs.^b) andU^b) 
are produced because the positive pulses hit the neigh- 
boring bubble when its internal pressure is not high. 



We have studied the interaction of microbubbles 
through pressure pulses to suggest a possible origin of 
the negative pulses found in broadband cavitation noise. 
The proposed scenario of noise generation is summarized 
as follows: When a large cavitation bubble collapses and 
fragments, a number of microbubbles are formed Tj. The 
microbubbles expand rapidly and then collapse to emit 
positive pressure pulses (the right half of Fig. [TO)) . The 
positive pulses hit and impulsively compress neighbor- 
ing bubbles to cause a brief but strong deceleration of 
the bubble surface. This deceleration, which creates a 
strong tension in the surrounding liquid, produces nega- 
tive reflected pulses (the left half of Fig. [TU)) , and then a 
signal consisting of positive and negative pulses is gener- 
ated. If the time-delay effect is considered, the positive 
and associated negative pulses are measured with a time 
interval determined by the relative position of the bub- 
bles and the sound speed or shock- wave velocity of water 
(which can highly increase in the vicinity of collapsing 
bubbles [2l|, [131). Spectral analysis has revealed that 
the frequency spectrum of the cavitation noise, particu- 
larly from the smaller bubble in a double-bubble system, 
becomes much more complicated and denser by bubble- 
bubble interaction. Though only a few negative pulses 
were observed in the double-bubble cases studied here, 
considering a larger number of bubbles and multiple scat- 
tering of sound may allow us to explain the large number 
of pressure pulses found in the recorded broadband noise. 

Discussing further details of the noise generation pro- 
cess, we have revealed that the amplitudes of the negative 
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reflected pulses depend on the instantaneous radius and 
the state of the internal gas of the reflecting bubble. In- 
terestingly, no negative pulses are generated when a sys- 
tem of identical bubbles is considered. This is because the 
positive pulse from a bubble hits the neighboring bubble 
just when it collapses, at which moment the bubble's in- 
ternal pressure is higher than the pressure of the incident 
positive pulse. This observation suggests that the surface 
of gas bubbles, a free surface, does not always produce 
a negative reflected pulse, and also that negative pulses 
cannot be described by theoretical models that consider 
only systems of identical bubbles. The presented find- 
ings would be useful to understand not only the shrimp 
emitted acoustic signal but also other types of cavita- 



tion noise found, e.g., in fluid machinery [3| and medical 
applications psl ]. 
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